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A derivative expansion of the effective average action beyond first order yields renormalization 
group functional flow equations which are used for the computation of critical exponents of the 
Ising universality class. The critical exponent v in D = 3 is consistent with high-precision methods. 



The effective average action 



Derivative expansion beyond first order 
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The effective average action T k [<p] ( see and USEE! ) 
is a functional of the fields (p. It interpolates between the 
classical action S at some microscopic ultraviolet scale A 
and the effective action T: T\ = S, lim^o Tfc = T- Its 
dependence on the momentum scale k is described by an 
exact renormalization group functional flow equation: 



d k T k = iTr{(lf 3 

,(2) 



Rk) 1 dkRk} 



where T k ' is the second functional derivative of Tk 
with respect to the fields and Rk denotes a momen- 
tum cutoff function with the properties lim^o Rk = 0, 
limfc^A Rk oo and lim g 2^ Rk(q 2 ) > ensuring the 
above limit properties of Tk . Universal properties of the 
renormalization group flow, however, are independent of 
the actual choice of Rk ■ 

(2) 

The flow of T k and higher n-point-functions can be de- 
rived from the flow of Tk by functional derivation, e.g. 



d k T k 2) = Tr{lf )2 (lf > 
-^r { lf>(lf) 

(3) 



Rk) 



- 3 d k R k } 
~ 2 dkRk\ 



where the flow of T, would involve T 



and T k and so 



on. In order to end up with a closed system of functional 
flow equations one needs to do a finite truncation of Tk- 

In the following we concentrate on a one-component real 
scalar field ip : M D — > M. There are the invariants 
p = \<p<p and pq 2 — ^d^ipd^ip featuring the O(l) symme- 
try of the Ising universality class. Here the lowest order 
truncation of the effective average action reads 

r fe = / d D x{u k (p) + \d^d^z k + o(d 2 )} 

in terms of the effective average potential Uk(p) and a 
field independent wave function renormalization Z k . The 
effective average propagator 

if V) + R k (q 2 ) = U'k(p) + 2pU' k '{p) + q 2 Z k + Rk{q 2 ) 
gives the flow of the effective average potential 

dkUk(p) = 5 / -^fn ui k{ p )+2 pufCp)+q'z k +R k {q') 
and dkZk can be determined from d q 2dkT k 2 \q 2 ) | g 2_ . 



In a second order derivative expansion of the effective 
average action there are contributions from up to four 
nonvanishing momenta. These are parametrized by the 
effective average potential and one respectively three lin- 
early independent wave function renormalizations in first 
respectively second order: 



T k = J d D x{U k (p) 
+ Z A ,k{p) 
+ Z B ,k(p) 

+ ZcAp) 
+ ZdAp) 



^d^tpd^tp 
■Q(pd, i d l *(pd v (pd v <p 



This truncation has been the starting point for a recent 
investigation in 6]. There an additional field expansion 
of the potential and the wave function renormalizations 
in powers of p is done. Here we keep the full field de- 
pendence but instead go only one step beyond first order 
and take only Uk, Z A ,k and Zn,k into account due to the 
following reasoning: 

(2) 

While Zam and Z B .k directly contribute to T k and thus 
enter all flow equations as q 2 and q 4 terms in the de- 
nominators of the integrands, Zb .k and Zc,k only give 
additive contributions to the four-point-vertex and the 
three-point-vertex respectively. Also, their flow requires 

d k ^k^ ana - ^fc^l . This is why we neglect Z B , k and Zc.k 
in this paper which corresponds to a truncation beyond 
first and below second order. 



We hope that this improves the critical exponent v com- 
pared to first order calculations. On the other hand, the 
anomalous dimension r] is related to the momentum de- 
pendence of the wave function renormalization and will 
probably suffer from the missing contributions. 



Substituting U := Uk, Zq := 
the truncation thus reads: 



Z A<k and Zi := ^ipipZ Dtk 



T k = Jd D x{U(p) 

+ Z (p) ■ IdpipWip 

+ z x {p) ■ ±d^<pd u d»<p + o(d 4 )} 

and the next step is to derive dkU, dkZo and d k Z\\ 
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Renormalization group flow equations 



Scale invariant flow equations 



Expanding the fields ip(x) = (p + x{ x ) m small fluctua- 
tions x{x) = ^2 e lqx x{q) in momentum space around 
a constant background (p yields the n-point-functions 



rW(gi,...,g n ) = 



Sf(qi) 6<p(q n ) 



r = 



T (2) (q, -q) = U' + 2pU" + q 2 -Z Q + q 4 -Z 1 

r( 2 ) (p, -p) = u' + 2 P u" + P 2 -z + P i -z 1 

rW(p,q,-p-q) = 

T (3) (p + «, -9, -p) = { 3(7" + 2pU"' 
+ (p 2 + {pq) + q 2 )-Z' 

+ (p 4 + 2p 2 (pq) + 3p 2 q 2 + 2(pq)q 2 + q 4 ) • Z[ } 
r( 4 ) (p, g, -g, -p) = W" + YlpU'" + 4p 2 U"" 
+ (p 2 + q 2 ) ■ {Z' Q + 2pZ'i) 
+ (p 4 + Ap 2 q 2 +q 4 )-(Z' l +2pZ'{) 

These are inserted into 

Q, r _ 1 f d D p d k R(p 2 ) 



(2tt) d r< 2 )(p,- P )+ii(p 2 ) 

d k TM(q,-q) 

— f r (3) (p,q,-p-^)r (3) (p+9,-9,-p) d h R{p 2 ) 

~ J (2k) d 



(2tt)° (rW(p2 ) + fl( p2 )) 2 (r (2) ((p+9) 2 )+i J ((p+9) 2 )) 

_ 1 f d D p T^(p,q-q-p) d k R(p 2 ) 
2 J (2tt) d (r( 2 ) (p 2 ) + R(p 2 )) 2 

and yield the flow equations 



8 k U = 8 k T 
dkZo 



In order to be able to perform the g 2 -derivatives one has 
to expand the denominator containing (p + q) 2 . Let N = 
r< 2 )(p 2 ) + R(p 2 ) and M = r< 2 )((p + q) 2 ) + R((p + q) 2 ): 



As we are interested in critical phenomena of second or- 
der phase transitions corresponding to fix points in the 
renormalization group flow, we write the flow equations 
in explicitely scale invariant form. 

To this end one substitutes zq(u) = Zo(p) / 'Zo(pk) ( where 
p k is the running minimum of U(p) ), a = Z (p k )k 2 ~ D p, 
u = k~ d U, t = — Lnfc, rj = —d t LnZ a (p k ) etc. such that 



d t u' = (-2 + ri)u' + (D-2 + r))au" 

- co • / 2 ° 

- ci • I\ 



d t z a = T}Z 
+ 4<jc c q 
+ 8erc ci 

+ 4<7CiCi 

+ 8cr c c 2 

+ 4(TClC 2 

+ 4ac 2 c 2 

- c 4 ■ / 2 ° 

- c 5 • 47 1 

St«i = (2 - 

+ 2(7C Co 

+ 4(TC Ci 
+ 2ac\C\ 
+ 4<jc c 2 

+ 4(TClC 2 
+ 2(7C 2 C 2 

-c 5 -7 2 ° 



( 



(D - 2 + r?)az 



/O+LD tO 

V D 3 
/1+2_D rl 

V D 3 
/0+3D rl 

V D 3 
/ 4+8T> t2 

V D 3 
/4+6£> t-3 

V S 3 



2+ D rl 

D ^4 1 D 
4+ _D 7 2 i 2 
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1 + f K P 



12+2D t3 
D 4 

J 4 + 



77)21 + (D - 2 + 77)0-^ 



0+ 2D j-0 
S J 3 ' 

0+ 2Z3 7-0 
5 J 3 

4+ 8D rl 
5 J 3 " 

16+22D 7-2 
z 3 



0+ 2D tO 1 
D J 4 + 
10+ 4D 7-1 
5 J 4 
8+ 6-D rl 



8+-D 
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4 + 



20+D 



^4°) 
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D 



D 

32+ 8D t2 
D J 4 " 

72+12D 7 3 
D J 4 



-l&D+D' 2 ts2\ 
D(D+2) 4> 

64+30D+-D 2 7^3\ 
D(n+2) Xi 4; 

96+42D+Z3 2 t^4\ 

TJJTJ+T) 4/ 



M TV 



( 2 pg+<7 2 )jV + (2 Pg + g 2 ) 2 AT 2 



TV2- 



Tv 3- 



"Tv 2- 



+ ... 



Then the mixed scalar products (pq) can also be elimi- 
nated: 

I / (gW, ? 2 ) • (pq) 2 = \ I ^f(p 2 ,q 2 ) p ^f 

\ I (Brf(p 2 , q 2 ) ■ (pq)* = \ J ^.f(p 2 ,q 2 )T^ 
\ I (£wf(P 2 > I 2 ) ■ (pq) k = for all odd k 
This finally gives the explicit flow equations: 
d k U = l-l[± 

- r " 4 p 2 jV 2 _ 2 p 2 jv ] , 

D N 5 D IV^J + ■•■ 
+ ... 



a fe z = 4pC 2 • 7 
a fc Zx - 2 P c 2 ■ 1 



N 



2N 2 
TV 5 " 



N 
TP 



f(P 2 )]-U£^f(P 2 )d k R(p 2 ) 



2 J (27T) 

C = 3J7" + 2pJ7'" 



where for abbreviation we have used the constants 

w = u' + 2au" 
co = w' 
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Cl 
C2 = Z\ 



C3 = w 
c 4 = z'n 



- 2<rw" 
2az'J 



c 5 = z[ + 2az'{ 
and the momentum integrals 



(p 2 



J™(w, z a ,z 1 ,r))[r] := 5 / (5^pT (P (L+^o+p^i+r)" ' 

d p p (p 2 )" 1 2(zo+2p 2 Zl +r) 2 gtr 



If d D p (p 2 ) m (2 Z1 +r) Btr 
2 J (2ts) u (w+p 2 z +p i z 1 +r)'-' 



which will be discussed next: 
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Momentum integrals 



Critical exponents 



The momentum integrals I, J and K can be computed 
efficie ntly in case of a linear cutoff. For the choice 
HESE3 Rk = Z k (k 2 - q 2 )S{k 2 - q 2 ) one has 



r = (1 -p 2 ) ■ 8(1 -p 2 ) 



r = -9(1 - p 2 ) - (1 - p 2 ) ■ 6(1 - p 2 ) 
r = 2 S(l - p 2 ) + (1 - p 2 ) ■ 8' [l - p 2 ) 
Btr = *§ffi = (2 - n + p 2 7 1 ) • 9(1 - p 2 ) 

and I, J and K are reduced to a single integral : 

C = (2 - fj W 

J™ = (z - 1)1? 
+ 2z!/™ +1 

X™ = 2(z - l) 2 /™ 1 

+ 8(zo-l)^V 1 1 

+ 8z 2 /" + i 2 

- v D (w + z + zi)~ n 
where H is a special case of a hypergeometric function 

H™(w,z ,zi) 

1 r d c p (p 2 )'" e(i-p 2 ) 

2 J (27r) T! («)+p 2 zo+p 4 2i+r)n 

1 

= v D J dx ((w+1)+{z *_ 1)x+ZlX 2 )n 







, fp fj- T a y z fc (n+fc-2-l)!fc! c ' 

(tu+1)" J Z^ fe! Z; (n-1)! 2! (fc-22)! 
fe=0 1=0 

oo fc/2 1 

«o (n+fc-2-1)! c' b k - 21 { j^^a+k 

(tu+l)" ^ ^ (n-1)! 2! (fe-22)! J u,i x 

fc=0 2=0 



fc/2 



t; D y^ y^ (n+fc-2-1)! c' fc fc - 2 ' 1 
(iu+1)™ Z^ Zv (n-1)! 2! (fe-22)! a+fc+1 



fe=0 2=0 



Starting with a quartic potential u = \(p — ka) 2 the 
flow equations are discretized on a grid and numerically 
solved for different values of ka- During the evolution 
towards k — > 0, the running minimum of the potential 
( u'(k) — ) may either end up in the massless sponta- 
neously broken phase ( m = 0, k — > oo ) or in the massive 
symmetric phase ( m 2 ~ k 2 u'(0), n = ). 

The critical k c leads to a second order phase transition 
characterized by fix points of all couplings and vanishing 
masses. Fine tuning ka around k c then yields the critical 
exponent u through the scaling law m ~ \T — T c \ v and 
the established proportionality |ka — k c \ ~ \T — T c \. The 
critical exponent r\ may be identified with the fix point 
■rf of the anomalous dimension. 

The critical exponents for the universality class of the 
three dimensional Ising model obtained in this way are 
compared in Table 1 to values obtained in leading order 
( Z\ = 0, Zo = 1, also known as 'local potential approx- 
imation' ), lowest order beyond leading order ( Z\ = 0, 
Zo = const, here called 'lowest order' ), and full first 
order beyond leading order ( Z\ = 0, here called 'first 
order' ). 

The exponents are in agreement with former calculations 
using the effective average action. For brevity we quote 
only results which are based on simulations as similar as 
possible. Still they differ in that they either use a field 
expansion ( an expansion of u and z in powers of p ) or a 
different infrared regulator or both. This might explain 
the numerical differences. 

For comparison we also give results from different ap- 
proaches. An extensive review of different techniques and 
recent results can be found in |l3| . We adopt their over- 
all estimates of exponents and errors. In addition we 
quote some of the newest and precisest results from high 
temperature expansions, perturbation series at fixed di- 
mension, e-expansions and Monte Carlo simulations plus 
some additional references. 



with 



-i = 2 d + 1 tt d / 2 T(D/2) 



a = m + y — 1 
b= i=a 

l + 'Uf 

C = 21- 

L 1+w 



The series expansion and reordering is valid and conver- 
gent as long as 1 — bx — cx 2 ^ for x 6 [0; 1]. Which is 
the the effective average propagator stays finite 

at all times. 



Conclusion and outlook 

The critical exponent v has significantly improved com- 
pared to first order. This is the first calculation using a 
linear cutoff that is consistent with high precision meth- 
ods. On the other hand, the anomalous dimension is 
sensitive to the neglected contributions as can be seen in 
comparison to 0. Its error has even increased in com- 
parison to first order. 

Hence the inclusion of the neglected wave function renor- 
malizations should be the next step. 
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V 


V 


effective average action to leading order ( this paper, Z\ = 0, Zq = 1 ) 


0.6491 





effective average action to lowest order ( this paper, Z\ = 0, Zq = const ) 


0.6271 


0.1120 


effective average action to first order ( this paper, Z\ — ) 


0.6255 


0.0503 


effective average action beyond first order ( this paper ) 


0.6303 


0.0793 


effective average action to leading order ( field expansion [11] ) 


0.6496 





effective average action to lowest order ( field expansion [12] ) 


0.6262 


0.1147 


effective average action to first order ( exponential cutoff [14] ) 


0.6307 


0.0470 


effective average action to second order ( exp. cutoff, field exp. [6] ) 


0.632 


0.033 


literature values from the review [13] 


0.6301(4) 


0.0364(5) 


25th-order high-temperature expansion [15], see also [16] 


0.63012(16) 


0.03639(15) 


seven- loop perturbation series at fixed D=3 [17], see also [18] 


0.6303(8) 


0.0335(6) 


five-loop order e expansion with boundary conditions 


0.6305(25) 


0.0365(50) 


Monte Carlo simulation [12], see also [20. 21] 


0.6297(5) 


0.0362(8) 



TABLE 1: Critical exponents of the three dimensional Ising universality class 
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